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Abstract 

Evolution of perturbed embedded solitons in the general Hamiltonian fifth-order Korteweg- 
de Vries (KdV) equation is studied. When an embedded soliton is perturbed, it sheds a one- 
directional continuous-wave radiation. It is shown that the radiation amplitude is not minimal 
in general. A dynamical equation for velocity of the perturbed embedded soliton is derived. 
This equation shows that a neutrally stable embedded soliton is in fact semi-stable. When the 
perturbation increases the momentum of the embedded soliton, the perturbed state approaches 
asymptotically the embedded soliton, while when the perturbation reduces the momentum of 
the embedded soliton, the perturbed state decays into radiation. Classes of initial conditions to 
induce soliton decay or persistence are also determined. Our analytical results are confirmed by 
direct numerical simulations of the fifth-order KdV equation. 

1 Introduction 

Embedded solitons are solitary-wave solutions of nonlinear evolution equations which reside at 
discrete points inside the continuous spectrum of the linear wave system. The existence of such 
waves has been known in various physical systems such as the fifth-order Korteweg-de Vries (KdV) 
equations [1, 2, 3], extended nonhnear Schrodinger equations [4, 5], coupled KdV equations [6], 
second-harmonic-generation (SHG) system [7], massive Thirring model [8, 9], three- wave system 
[10], and many others [11, 12]. In [7], such waves were given the name "embedded solitons", and 
their distinct semi-stabihty property was revealed on heuristic ground. This semi-stability means 
that when a perturbation increases a certain positive-definite quantity (energy or momentum) 
associated with the embedded soliton, then the perturbed state approaches asymptotically the 
embedded soliton. However, when the perturbation decreases energy (momentum) of the embedded 
soliton, the perturbed state decays into radiation. 

The semi-stability property was later proved rigorously for embedded solitons in the perturbed 
integrable fifth-order KdV equation [13] and in the generalized SHG system [14]. The method in 
[13] follows the soliton perturbation technique and describes embedded solitons as critical points of 
a first-order dynamical system. Both location and stability of critical points can be studied within 
the reduced system. The other paper [14] develops the normal form analysis which relies on the 
known existence and linearized stability properties of embedded solitons. It proves the nonlinear 
semi-stability of embedded solitons through wave resonance mechanisms. The latter method does 
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not rely on the integrability of the original system, and it can be extended to any embedded-soliton- 
bearing system under certain assumptions. 

Nonlinear semistability of embedded solitons is an interesting phenomenon because it occurs 
beyond the linear stability. The linear stability of solitary waves in the fifth-order KdV equation was 

studied by using the energy-momentum methods [15, 16] and the symplcctic Evans matrix methods 
[17]. A solitary wave is linearly unstable if the linearization operator possesses an eigenvalue with a 
positive real part. If all eigenvalues lie on the imaginary axis, the wave is called neutrally stable. A 
neutrally stable wave can however still be unstable due to algebraic instabilities [18]. For embedded 
solitons, the situation is different. Single- hump embedded solitons are generally neutrally stable in 
the linearized problem. However, a discrete zero eigenvalue of the linearization operator is embedded 
inside the continuous spectrum of this operator. Because of this, a nonlinear resonance between the 
embedded zero eigenvalue and the continuous spectrum can tunnel energy of a perturbed embedded 
soliton into continuous-wave radiation. The energy loss does not always destroy the embedded 
soliton though. If the perturbation increases the energy (momentum) of the embedded soliton, 
then the radiation becomes weaker and weaker as the perturbed state asymptotically approaches 
the embedded soliton. But if the perturbation decreases the energy (momentum) of the embedded 
soliton, the radiation becomes stronger and stronger, and the embedded soliton is then destroyed. 
Thus, the semi-stability of embedded solitons is an intrinsically nonlinear phenomenon beyond 
linear stability. 

Besides single-humped embedded soliton in the fifth-order KdV equation, multi-humped embed- 
ded solitons may also exist [3, 19]. However, they are typically linearly unstable [11, 12]. The 
nonlinear semi-stability may occur in the system only if the linear instability is suppressed. Thus, 
it makes sense to consider here only neutrally stable single-humped embedded solitons. 

We emphasize that semi-stability is not the same as instability. By controlling the energy (mo- 
mentum) of the initial perturbation, we can induce either asymptotic persistence of the embedded 
soliton or its rapid disappearance. This is an ideal mechanism for switching (quantization) appli- 
cations. 

In this article, we study the evolution of perturbed neutrally stable embedded solitons in a 
general Hamiltonian fifth-order KdV equation. We use the normal form analysis developed in [14] 
but simplify many statements and proofs. Wc show that when an embedded soliton in the fifth- 
order KdV equation is perturbed, it sheds the continuous-wave radiation in front of the soliton. 
The radiation amplitude is not minimal in general. We also derive the velocity equation for the 
perturbed state which proves the semi-stability property of embedded solitons. In addition, we 
determine what initial condition leads to soliton decay, and what initial condition leads to soliton 
persistence. Numerical simulations of the fifth-order KdV equation show excellent agreement with 
the analytical predictions. 

2 Tail amplitudes of symmetric nonlocal waves 

We consider the general Hamiltonian fifth-order KdV equation 

Ut + Uxxx+Uxxxxx+Wiu)]x = 0, (2.1) 

where the nonlinear term N(u) is of the form 

N{u) = aou^ + aiuuxx + <^2ul + a^u^. (2.2) 
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The equations above are Hamiltonian if and only if ai = 2a2 [13]. The conserved Hamiltonian 
functional H(u) is 



H{u) = / 

J — ( 



"0 3 1 2 ,1 2 "1 2 "3 4 

3 2^ 2^^ 2 ^ 4 



dx. (2.3) 



When Eq. (2.1) is Hamiltonian, it also conserves a positive-definite quantity 

/oo 
u^dx, (2.4) 
-oo 

which is interpreted as momentum. Whether Eq. (2.1) is Hamiltonian or not, it always conserves 
the quantity M(u) = J^^udx, which is interpreted as mass. 

We look for moving stationary solutions in Eq. (2.1) of the form 

u{x, t) = U{x - ct) = U{C), (2.5) 

where c is the wave velocity. Substituting Eq. (2.5) into (2.1) and integrating once, we obtain the 
ordinary differential equation (ODE) for U (^) as 

%«« + %-cC/ + iV(?7) = 0, (2.6) 

subject to zero boundary conditions at infinity: lim|^|^oo U = 0. When [/ ^ 0, Eq. (2.6) 
becomes linear, and its characteristic equation is quartic. The roots of this characteristic equation 
are :^ik and ibK, where 



. = *(.) = ^^l±i^, . = „(.) = ^^l±i^. (2.7) 

When c > 0, roots ibf/c are imaginary and ibK real; when > c > —1/4, all four roots are imaginary; 
when c < —1/4, all four roots are complex. The embedded soliton, if it exists, arises from a saddle- 
node bifurcation [19], when the stable-unstable manifolds (real roots) correspond to the exponential 
tails of the embedded soliton, while the center manifold (imaginary roots) correspond to the tails 
of the continuous-wave radiation whose amplitude vanishes. Therefore, the embedded soliton may 
exist in Eq. (2.6) only for c > 0. We notice that non-embedded solitary waves with oscillatory and 
decaying tails may exist in Eq. (2.6) for c < —1/4 [20] but such solutions are beyond the scope of 
this paper. 

Only symmetric embedded solitons of Eq. (2.6) are considered here. Based on energy flux 
consideration, it is generally believed that asymmetric embedded solitons of Eq. (2.6) do not exist 
[21, 22]. If the system (2.1) is the perturbed integrable fifth-order KdV equation, the non-existence 
of asymmetric embedded solitons in Eq. (2.6) was proved in [13]. 

When c > 0, the tail amplitude of symmetric nonlocal waves of Eq. (2.6) is given asymptotically 

as 

U{^-c,5)^r{c,6)s\n{k\^\+5), |C| ^ oo. (2.8) 

Here r is the amplitude of the oscillatory tail, and 5 is the tail phase. Embedded solitons can be 
found numerically by a shooting method for Eq. (2.6) with t/'(0) = U"'{0) = 0. Parameters of 
the shooting method U{0) and U"{0) satisfy only one condition that removes the exponentially 
growing term ~ e'^'^L It implies that one parameter in symmetric nonlocal wave solutions of Eq. 
(2.6) is free in addition to the wave velocity c. We choose this free parameter to be the tail phase 
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6. The tail amplitude r{c,6) may vanish at certain discrete velocity values c = ces (in such cases, 
the phase S becomes irrelevant). When it happens, we get embedded solitons. The sech^ embedded 
solitons in Eq. (2.1) has been investigated comprehensively in [1, 13], and a family of one, two, 
and oo embedded solitons have been shown to exist in different parameter regions. Single-humped 
and multi-humped solutions of (2.6) were studied also in [19]. We assume here the generic case of 
a co-dimension one bifurcation, when the embedded-soliton velocity ces is a simple zero of the tail 
amplitude r{c,6), i.e., 

r(c, S) = R{d){c - ces) + 0(c - ces)', (2.9) 

where 

R{S) = ^^{CES,S) ^0 (2.10) 

is the slope of the tail amplitude r of the symmetric nonlocal waves at the embedded-soliton velocity 
Ces S'lid phase 5. 

In the rest of this section, wc derive the analytical expression for the tail amplitude r(c, S) when 
the nonlocal wave velocity c is close to the embedded-soliton velocity ces- For this purpose, we 
expand the nonlocal solution ?7(^; c, 5) as a perturbation series: 

U{(;c,d) = UEs{0 + {c-CEs)Ui{^;5) + O{c-CEs)^. (2.11) 

When this expansion is substituted into Eq. (2.6), the function Ui{^;6) is found to satisfy the 
following inhomogeneous linear equation: 

CUi{C;S) = Ues{0, (2.12) 
where C is the linearization operator of Eq. (2.6) at embedded-soliton velocity ces, i-e.. 



= ^ + ^ - CES + 2aoL/ES + 3a3t/Es + «! ^ V d| 

Here the Hamiltonian condition ai = 2a.2 has been utilized. Note that the operator L is self-adjoint 
in the Hamiltonian case. 

In order to solve the inhomogeneous equation (2.12), we need to know homogeneous solutions. 
Operator L has four homogeneous solutions V'n (1 < < 4). The first solution is symmetric and 
bounded with asymptotic behavior: 

V^i(0 sin(fe,|^| + (^,), \i\ ^ 00, (2.14) 

where ky = A;(ces) is the resonant wavenumber and is the tail phase. The second solution is 
anti-symmetric and bounded with asymptotic behavior: 

^-2(6 sin(A:,e ± 5a), e ^ ±00, (2.15) 

where ba is the tail phase. The third solution is anti-symmetric and localized, ^^3(0 = ^es(0- is 
related to spatial translational invariance of the system (2.1). The last solution ?/'4(^) is symmetric 
and unbounded. 

Now we can solve the inhomogeneous equation (2.12) for \J\. Since we only consider symmetric 
nonlocal waves, the correction term f/i(^) is symmetric as well. Note that c?C/(^; ces, <5)/9c is an 
inhomogeneous solution of Eq. (2.12) for any 8 values. This can be seen by differentiating Eq. 
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(2.6) with respect to c, and then setting c as ces- Thus the general bounded symmetric solution of 
Eq. (2.12) can be written as 

UiiO = ^(C;cES,<5o) + 7iV'i(0- (2-16) 

Here Sq is any fixed phase, and 71 is an arbitrary constant. The homogeneous solutions ^^2(0 
ip3{^) are not included as they are anti-symmetric, while the solution t/j4:{0 unbounded. 

The asymptotic oscillatory behavior of a general Ui solution (2.16) can be obtained from Eqs. 
(2.8) and (2.14) as 

Ui{C;6)-^R{So)sm{kr\^\ + So)+iiMkr\^\+Ss), \^\ ^ 00. (2.17) 

where function R{S) is defined by Eq. (2.10). On the other hand, expanding the tail asymptotics 
(2.8) of the nonlocal wave U{^;c, 6) into a power series of (c — ces) and then comparing it with Eq. 
(2.11), we conclude that the asymptotic behavior of a general Ui solution should also be 

UiiO R{d)sm{kr\^\+d), ^ 00. (2.18) 

For convenience, we fix here 60 as Sq = Sg + 7r/2 and then equate the two asymptotics (2.17) and 
(2.18). As a result, we find that 

sm{6 — Os) 

and 71 = R{Ss + §) cot(5 — Sg). Consequently, to leading order in (c — ces), the tail amplitude r of 
symmetric nonlocal waves is 

r(c, 6) = :?%±|l(c - ces) + 0(c - ces)'. (2.20) 
sm{d — Os) 

From this formula, we conclude that at a given velocity c close to ces, the tail amplitude |r| of 
symmetric nonlocal waves is minimal at the value S = 6^1^ where 

(5s + mod(7r). (2.21) 

Formulae (2.20) and (2.21) are the main results of this section. 



3 Dynamics of embedded solitons under perturbations 

In this section, we study dynamics of a linearly neutrally stable embedded soliton under small 
perturbations. We assume that the linearization operator has no unstable eigenvalues. This as- 
sumption is necessary as the weak nonlinear semi-stability would be ineffective in the presence 
of strong linear instability. We also assume that the linearization operator has no discrete non- 
zero embedded eigenvalues. This assumption is necessary as the non-zero embedded eigenmodes 
may be in resonance with the continuous spectrum through nonlinear coupling, thus affecting the 
dynamics of embedded solitons. Thirdly, we assume that the zero eigenvalue corresponds to the 
localized eigenfunction C/es(0 ^^'^ algebraic multiplicity two with an associated eigenfunction 
dU/dc{^;S). For embedded solitons, the zero eigenvalue is always embedded into the continuous 
spectrum of the linearization operator. Eigenvalues that correspond to localized eigenfunctions 
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occur as zeros of the Evans function (a determinant of scattering coefficients) [17]. The algebraic 
multiphcity of eigenvalues is defined as the multiplicity of zeros of the Evans function. Thus, our 
last assumption is that the Evans function has a double zero at A = of the linearized operator. 

We will use below the internal perturbation analysis described in [14] (see also [18]). The idea 
is to recognize that under small perturbations, the cigenfunctions U'^g{(,) and dU/dc{^;5) for the 
double embedded eigenvalue A = of the linearized problem renormalize the location and velocity 
of the embedded soliton. For small perturbations, the velocity c{t) of the embedded soliton changes 
on a slow time scale. We will derive a dynamical equation for c{t) by separating the slow and fast 
changes in evolution of a perturbed embedded soliton. 

In the moving coordinate, 

^ = x— cdt — xo, (3-1) 
Jo 

the fifth-order KdV equation (2.1) can be written as 

ut - cu^ + u^^^ + u^^^^^ + [N{u)]^ = 0. (3.2) 

We expand the perturbed embedded soliton and its slowly varying velocity into the following 
perturbation series: 

n{^,t) = UEs{0 + MT)ui{^,t)+e^U2{^,t,T) + 0{e^), (3.3) 

and 

c{T) = CES + eci(r) + 0(6^), (3.4) 

where T = et, and e is a small parameter. At order e, we obtain the governing equation for «i(^, t) 
as 

uu + [/:ui]^ = U^siO, (3-5) 

where C is the same linearization operator as defined in Eq. (2.13). The initial condition for Eq. 
(3.5) can be obtained from Eq. (3.3) as 

wi(4,0 = — . 3.6) 

eci(O) 

The initial value for the soliton velocity ci(0) can be found by projecting the initial deviation 
tt(^, 0) — f7Es(0 onto 9C/(^; ces, (5a)/9c, where the phase 5^ is given by (2.15). The projection is 
based on the spectral decomposition developed in Appendix A (see (A. 4)) and is given by 



roof^Es(e)f (C;CES,5a)c/C 



eci(o) = Too ;; a;;^,: — ""i^- (3-7) 



Next, we solve the inhomogeneous equation (3.5). We adopt a less formal but more intuitive 
approach here. A more rigorous calculation of the same results is presented in Appendix A. 

The inhomogeneous term in Eq. (3.5) acts as a driving localized force. The homogeneous part 
at large j^j values supports oscillatory solutions with wavenumber kr- Thus due to forcing and 
resonance, these oscillatory tails will be excited over time. The group velocity of these oscillatory 
tails in the moving frame (3.1) can be found from the dispersion relation as Cgr = 2k^{2k^ — 1), 
which is always positive since ces > and kr > 1. Thus, these oscillatory tails always appear ahead 
of the embedded soliton. Behind the embedded soliton, there is the possibility that a fiat shelf may 
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develop as in the perturbed KdV equation [23, 24] (see also [18]). If a shelf develops, it moves to 
the region x <C — 1 at the velocity — ces in the moving coordinate system (3.1). 

Thus, at large times t » 1, the boundary conditions for the solution ui{$,,t) are 

ui{^, t) I ^-^f "(^'-^ + ^-d)i^(cg.i -0, e » 1, (3.8) 
^ \ RoH{-C-CEst), ?<-l, ^ ^ 

where i?rad is the oscillatory-tail amplitude, (^jad is its phase, Rq is the height of the trailing shelf, 
and H{x) is the step function, i.e., H = 1 when a; > 0, and H = otherwise. Below, we determine 
the tail amplitude i?rad and the phase Jrad- We also show that the shelf is not excited in the present 
situation, i.e., Rq = 0. 

Our calculations of Rq, i?rad and 5rad are based on the observation that, as t goes to infinity, the 
transient part of the solution ui{^,t) for Eq. (3.5) dies out, and ui{^,t) approaches a steady state 
uis{$,) where 

This steady-state solution satisfies the same equation (3.5) except that the time derivative in (3.5) 
is dropped, i.e., 

[ji:uu]^ = u^s{0- (3.10) 

Integration of this equation with respect to ^ gives 

jO-uu = Ues + V, (3.11) 

where ?7 is a constant. To determine rj, we substitute the boundary condition (3.9) of solution 
uisiO at ^ » 1 into Eq. (3.11) and find that r/ = 0. Then substitution of the boundary condition 
(3.9) at ^ <C — 1 into Eq. (3.11) readily shows that Rq = 0. Thus, the fiat shelf is not excited in 
the present situation. 

Since Ro = r) = 0, the inhomogeneous equation (3.11) for uis{$,) becomes the same as Eq. (2.12). 
The general bounded solution for uis(^) is 

uisiO = ^(?; CES, <^o) + riVi(0 + ^2M0, (3.12) 

where Sq is any fixed phase and ri^2 are constants. The homogeneous solution V'3(0 is excluded 
by a simple position normalization of the embedded soliton. The boundary condition of solution 
(3.12) at infinity can be obtained from Eqs. (2.8), (2.14) and (2.15) as 

(c\ / R{6o)sm{kr^ + 5o) +Tism{kr(, + 5s) + T2sm{kr^ + 5a), ^ oo, . . 

This boundary condition should match condition (3.9) with Rq = 0. For convenience, we fix So = da- 
Then the matching condition gives the radiation amplitude i?rad and phase t^rad as 

iJrad = 2R{5a), ^rad = 5a, (3.14) 

and 

Ti = 0, r2 = R{5a). (3.15) 
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Formulae (3.14) are important results of this section. They show that the radiation phase 5rad is 
equal to the phase Sa of the anti-symmetric homogeneous solution ■02(0 [s^e Eq. (2.15)], while the 
radiation amplitude Rrad = 2-R((5a). Since the minimal tail amplitude of symmetric nonlocal waves 
occurs at phase 5„ii^ = [5s + §) mod(7r) [see Eq. (2.21)], and 6a ^ (Ss + f) mod(7r) in general (see 
Sec. 4 for an example), we conclude that the radiation amplitude Rj-ad generally is not minimal, 
i.e. Rj-ad 7^ 2i?((5min)- In the numerical work for the KdV equation plus the fifth-order derivative 
[25] (Sec. 16.6), it was mentioned without proof that radiation tail amplitude was minimal. That 
statement does not agree with our general analysis. But if the fifth-order KdV equation is integrable, 
then the relation da = 6s + ^ holds [13], i.e. the radiation amplitude is indeed minimal to the leading 
order of the perturbation theory for nearly integrable fifth-order KdV equations. 

When the first-order solutions (3.12) and (3.15) are substituted into the perturbation expansion 
(3.3), the solution can be re- written as 

u{C, t) = {U{C; c, 5a) + (c - CEs)R{Sa)MO + 0[(c - CEsf]}H{c^t - 0, t » 1, (3.16) 

where c{T) is given by Eq. (3.4). This solution describes the slow evolution of the perturbed 

embedded soliton in the fifth-order KdV equation (2.1), while the fast radiation part produced by 
a general initial condition for ii(^,0) is neglected in the asymptotic limit i 3> 1. Solution (3.16) 
up to order 0(c — ces) consists of a symmetric nonlocal wave U{^;c,5s,) and an anti-symmetric 
term "02 (0- T^^^s anti-symmetric term is generated in the initial- value evolution problem due to 
the (radiation) boundary condition (3.8) with Rq = 0, i?rad = 2i?((5a), and 5j-ad = ^a- Since function 
ip2{C) is anti-symmetric, the radiation amplitude is canceled behind the embedded soliton and is 
doubled ahead of the soliton. It is also noted that ^^2(0) = 0, thus the amplitude of solution (3.16) 
at soliton center ^ = is then the same as that of the symmetric nonlocal wave i7(^;c, (5a). This 
fact will be used in Sec. 4 in our comparison between the analytical and numerical results on the 
amplitudes of perturbed embedded solitons. 

When the radiation amplitude -Rrad and phase J^ad ^ire found, we are ready to derive the dynamical 
equation for the velocity c(T) of a perturbed embedded soliton. This equation can be derived 
in several different ways (see [13]). The simplest way is to use the local or global momentum 
conservation law when the system (2.1) is Hamiltonian. The derivation using the global momentum 
conservation law (2.4) is presented below. The derivation using the local momentum conservation 
is contained in Appendix B. 

To derive the velocity equation, we substitute the perturbation expansion (3.3) into the momen- 
tum integral (2.4). When terms up to order are retained, we get 

-j [ul^ + 2eUEs{ciUi + eu2) + e^clul}di = Q. (3.17) 

Keep in mind that solutions ui{^,t) and U2i6.,t) at the center region ^ ~ 0{1) become stationary 
as t S> 1. As a result, the term involving U2 in Eq. (3.17) can be dropped because the integral of 
the product f/Es(0^2(C)0 becomes constant at large times. The stationary solution uis{^) is given 
by Eqs. (3.12) and (3.15). Thus, 

roc roc QJJ 

E= UEsuisd^= i7Es(0^(e;cES,W- (3-18) 

J ~oo J ~oo (JC 

Lastly, the solution for ni(^,t) develops an oscillatory tail ahead of the embedded soliton. This tail 
has amplitude i?rad given by Eq. (3.14), and it moves at its group velocity Cgr. The average energy 
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(iif) of the sinusoidal tail is i-^-' 2-R^((^a)- Thus, 

V=-j_^ ujd^ = 2R\5a)cgr = ^kl{2kl - l)R\5a). (3.19) 

When relations (3.18) and (3.19) are substituted into the momentum equation (3.17) and eci 
replaced by c — ces [see Eq. (3.4)], we finally obtain the dynamical equation for the perturbed 
embedded-soliton's velocity c as 

- = -/3(c-CEs)^ (3.20) 
where the coefficient (5 = T /2E. The solution of Eq. (3.20) is 

where cq is the initial condition for velocity c{t). The formula for cq can be obtained from Eqs. 
(3.4) and (3.7) as 

r^^Es(OKe.o)-^Es(C)]c ^C 
root^Es(6f (C;CES,W 
Once the initial perturbed embedded soliton "^(^,0) is specified, then cq is fixed as above. 

The asymptotic equation (3.20) is the key result of this paper. When /? > 0, this equation 
shows that its fixed point c = ces is semi-stable: any perturbation with cq > ces is stable, and 
any perturbation with cq < ces is unstable. When translated into the original partial differential 
equation (2.1), it means that the embedded soliton is semi-stable. Depending on the type of initial 
perturbations, the embedded soliton can persist, or be destroyed. 



Finally, when E vanishes, the zero embedded eigenvalue has multiplicity higher than two, which 
results in linearized (algebraic) instability of the embedded soliton (see, e.g., [18]). We have excluded 
such linearized instability in our assumptions above. 



4 Comparison with direct numerical simulations 

In this section, we directly simulate the original partial differential equation (2.1), and compare the 
results with our analytical theory above. The system parameters we choose are 

ao = 5, ai = 5, a2 = 2.5, 0:3 = (4-1) 

in Eq. (2.2). Note that these parameter values are equivalent to ao = ai = 1, a2 = 0.5 and as = 
after variable u and time t are rescaled. At these parameter values, Eq. (2.1) is Hamiltonian. The 
fifth-order long-wave model equation studied by Champlneys and Groves [3] corresponds to our 
equation (2.1) with ao = 1 and 0:3 = 0. We have also tested other parameter values with ai = 2a2 
and found similar results. For instance, in the third-order Hamiltonian long-wave approximation to 
the water-wave problem as derived by Craig and Groves [26], the parameter values (after variable 
rescaling) are ao = 1, ai = ~|) 0:2 = — | and as = 0. Comparison between our theory and numerics 
for this set of parameters is qualitatively the same as that for the parameters (4.1). 

With the parameters (4.1), the fifth-order KdV equation (2.1) has an embedded soliton 

Ues {x,t) = 0.9 sech^ f VOS (x - cesOI (4-2) 



9 



at the exact wave speed ces = 2.64 (see [11]). The approximate phase values Sg and 6a in the hnear 
modes V'1,2 of the linearization operator C are found numerically (by the shooting method) as 

6s = 2.1815, 6a = 0.5737. (4.3) 

Note that the difference between these two phases here is not equal to n/2, thus the radiation tail 
amplitude in perturbed embedded solitons is not minimal. However, this phase difference differs 
from 7r/2 only by 0.037. Thus the radiation tail amplitude is rather close to its minimal value. 
At phase 6 = 6a (which is the radiation phase), we have numerically obtained the tail amplitude 
curve r(c, 6a) of symmetric nonlocal waves from Eq. (2.6) at various velocity c values, again by the 
shooting method. The results are shown in Fig. 1. As expected, the tail amplitude r is non-zero 
for c > except when c = ces- The slope R{Sa) at embedded-soliton velocity ces is found to 
be R{6a) = —0.0652. The kr value can be quickly obtained from Eq. (2.7), and the integral in 
Eq. (3.18) can be readily determined numerically. Prom all these values, we finally found that 
/? = 0.0868. With this /? value, our analytical formula for the velocity of perturbed embedded 
solitons is then given by Eq. (3.21). 

In order to verify our analytical theory, we have numerically simulated the original wave equation 
(2.1) with system parameters (4.1) and initial condition 

uix,0)=hUEsix,0), (4.4) 

where h is a constant coefficient. Note that h = 1 gives the exact embedded soliton, and h ^ 1 
gives a perturbed embedded soliton. Our numerical scheme is the integrating factor method as 
described in [27]. The x interval is taken as 400 units long, and 1024 grid points are used. The 
time stepsize is 10~^. To prevent radiation from re-entering the simulation region through periodic 
boundary conditions, we have used a damping condition near the boundaries. In our simulation, 
we have also adopted a frame moving at the cmbcddcd-soliton's velocity ces (but the results will 
still be presented in the original frame). Our numerical scheme has been tested with the exact 
embedded soliton (4.2) as the initial condition. It has also been tested on a related system — the 
integrable fifth-order KdV hierarchy equation. Furthermore, we have tried different grid points and 
time stepsizes. These tests show that the numerical error in our scheme is on the order of 10~^. 

We have run our numerical scheme on two typical initial conditions (4.4) with h = 1.05 and 0.95. 
The results are presented in Figs. 2 and 3 respectively. In the former case, the perturbed state has 
momentum P higher than the embedded soliton's [see Eq. (2.4)]. Because of this, the perturbed 
state initially moves a little faster than the unperturbed embedded soliton [see Fig. 2(b)]. But its 
speed as well as amplitude slowly decrease due to continuous wave radiation which moves ahead 
of the main pulse [see Fig. 2(b,c)]. This tail radiation at t = 20 is shown in Fig. 2(a) (note 
that the tail decay near the right end of the x-interval is due to our damping boundary condition. 
The actual tail length is much longer). But the tail amplitude decreases also in the process [see 
Fig. 2(d)]. Thus energy radiation is decreasing. Eventually the perturbed state asymptotically 
approaches the unperturbed embedded soliton, which is clearly seen in Fig. 2(b,c). 

When h = 0.95, the perturbed state has momentum P lower than the embedded soliton's. In 
this case, due to continuous wave radiation which intensifies over time [see Fig. 3(a,d)], the speed 
and amplitude of the perturbed embedded soliton both decrease well below their corresponding 
values of the unperturbed embedded soliton [see Fig. 3(b,c)]. When the amplitude of the main 
pulse has dropped significantly, it can no longer sustain high radiation tails. Thus tail amplitudes 
start to decrease [see Fig. 3(d)]. Eventually, the embedded soliton is destroyed by perturbations. 
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The above numerical simulation results agree both qualitatively and quantitatively with our 
analytical theory. Qualitatively, when h = 1.05, as the initial velocity is above ces, formula (3.21) 
predicts that the pulse velocity will asymptotically approach ces; when h = 0.95, the velocity will 
decay far below ces- This semi-stability behavior is accurately reflected in the numerical results. 
Quantitatively, we have also compared the pulse velocity, amplitude and tail amplitude of analytical 
predictions to those of numerical results. The analytical prediction for pulse velocity is given by 
formula (3.21). The initial condition cq is calculated from formula (3.22). We found that when 
h = 1.05, Cq ~ 2.775, and when h = 0.95, cq ~ 2.515. The analytical prediction for the pulse 
amplitude is the center amplitude of symmetric nonlocal waves at analytically predicted velocity c 
[see Eq. (3.16)]. The analytical prediction for tail amplitude is (c — CEs)-^?radj where Rrad is given 
by Eq. (3.14). These analytical predictions have been plotted in Figs. 2(b,c,d) and 3(b,c,d) as well 
for comparison. In the case h = 1.05, the quantitative agreement between theory and numerics 
is excellent at all times. In the other case h = 0.95, the quantitative agreement is good at the 
beginning, and gets worse at larger times. The good agreement in the former case is because the 
main pulse remains close to the embedded soliton at all times, thus the perturbation theory works 
well. In the latter case, the main pulse deviates significantly from the embedded soliton at large 
times. When that happens, the perturbation theory breaks down. 



5 Conclusion 

In this article, we have studied the evolution of perturbed embedded solitons in a general Hamilto- 
nian fifth-order KdV equation (2.1). We have shown that when an embedded soliton is perturbed, 
it sheds continuous-wave radiation in front of the embedded soliton. The amplitude of this con- 
tinuous wave is not minimal in general. Behind the embedded soliton, no flat shelf is created. We 
have further derived the velocity equation of a perturbed embedded soliton. As a result, the semi- 
stability property of embedded solitons is analytically proved. In addition, we have obtained the 
conditions under which a perturbed embedded soliton will decay or persist. We have also simulated 
the fifth-order KdV equation numerically. The numerical results agree well with the analysis both 
qualitatively and quantitatively. 

The analysis and the final dynamical equation (3.20) are similar to those found in [14] for gener- 
alized SHG models. Thus, in spite of differences in the spectral properties of linearization operators 
in these two models, the nonlinear resonance between the embedded soliton and the continuous- 
wave radiation has common features under assumptions listed in Section 3. Obviously, the same 
method can be applied to any other embedded-soliton-bearing Hamiltonian system. In fact, the 
system does not even have to be Hamiltonian. A non-trivial conservation law such as power or 
momentum would be sufficient to guarantee the semi-stability property of embedded solitons (see 
[14]). The open problems going beyond the present study include generation and collisions of sev- 
eral embedded solitons, as well as further engineering applications of embedded solitons in applied 
science. 
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Appendix A 



In this appendix, we present a more rigorous approach for solving the inhomogeneous equation 
(3.5) for the first-order solution t) by using a spectral decomposition method for the linearized 
problem, 

{Ccl){^;k))i: = in{k)(l){e,k), (A.l) 

where Q,{k) = k{k'^ — k^ — ces) , and (p{^; k) are continuous- wave eigenfunctions normalized according 
to the boundary condition: 

^(^;fc)^e^^^ as (A.2) 

The potential terms with C/es(C) decay exponentially at large |^| in the operator C given by Eq. 
(2.13). The eigenfunctions (p{(^,k) may have up to three Fourier oscillatory terms in the limit 
^ +00, which match with the roots of the equation: Q{k) = Q. We will compute asymptotically 
the Fourier-type integrals [see Eq. (A. 10) below] at the resonant values k = ±kr, where kr = k{cEs) 
and k{c) is given by Eq. (2.7). At the resonant values, the eigenvalue parameter Q is zero, i.e. 
^}{±kr) = 0. It could be found from Eqs. (2.14) and (2.15) that the boundary condition for 
(/>(^, zizkr) in the limit ^ — > -|-oo is: 

(f>{^;±kr)^a±e^'''^^ + b±e'^'''^^ + c±, as ^ » 1, (A.3) 

where 

sm{ds-da) sm[ds - da) 

Under assumptions described in Section 3, the solution for ui{$,,t) can be decomposed through 
eigenfunctions of the linearized problem (A.l): 

ni(e, t) = — (C; CES, <5) + wik, t)0(C, k)dk + aU^siO, (A.4) 

where a is constant. The first term in (A.4) solves the inhomogeneous part of Eq. (3.5). Since 
the double eigenvalue O = is embedded into the integral at the resonant points k = ±kr, the 
inhomogeneous term is not independent and can be decomposed through the same eigenfunctions: 



— (C;cES,(^) = / .2 1,2 dk, (A.5) 



where F{±kr) 7^ 0. The singular (pole) part in (A.5) describes the non-localized oscillatory tail 
(2.18) as 1^1 » 1. With Eq. (A.5) substituted into Eq. (A.4), at time t = 0, Eq. (A.4) becomes 

/oo 

w{k,0)cl)i^,k)dk + aU^siO, (A.6) 
-00 

where 

w{k,0) = -^^^+wik,0). (A.7) 

When the initial condition mi(^, 0) is localized, then the spectral component u;(fc,0) is free of pole 
singularities. 

The complex amplitude w{k,t) in Eq. (A.4) satisfies the trivial evolution equation: 

-^+in{k)w = 0. (A.8) 
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Solving the initial-value inhomogcncous problem (3.5) with the spectral decompositions (A. 4) and 
(A. 5), we find the following integral representation for ui{^,t): 

M^,t) = r (i - e-^W*) dA; + r u)(fc,0)e-W^(^,fc)dfc + «c/^g(^). (A.9) 

J —oo fC fvj, ^ ' J —oo 

The second term in Eq. (A.9) represents the non-singular part produced by the initial condition 
ui{^, 0). The first term in Eq. (A.9) represents the singular (pole) part produced by slow evolution 
of the embedded soliton. The singular integral term describes the nonlinear resonance between the 
embedded soliton and the continuous-wave radiation. 

The singular (pole) term in Eq. (A.9) occurs for k = zt/cj., i.e. when i}{zizkr) = 0. This term 
represents the oscillatory-tail radiation that diverges from the embedded soliton with the group 
velocity Cgr = ^l'{kr) = 2k^{2k^ — 1). Since ces > and k^ > 1, then Cgr > 0, i.e. the oscillatory- 
tail radiation occurs ahead but not behind of the embedded soliton. We prove this conjecture by 
using the pole decomposition technique [14]. The singular contribution from the integral (A.9) can 
be evaluated in the asymptotic region ^ <C — 1, t ^ 1 such that is constant: 



[F{kr)e""'^ - F{kr)e-'^^^) [sign(C/t) - sign(C/t - c^)] = 0. (A.IO) 



Here we have used the boundary condition (A. 2) and the symmetry relation F{kr) = F{—kr). 
Similar but lengthy computations of the integral (A.9) with the boundary condition (A.3) in the 
region ^ » 1 prove that the boundary conditions for the solution ui{^,t) of Eq. (3.5) in the limit 
t S> -|-oo and ^/t constant is 

„ .^ i -Rrad Sm{kr^ + (5rad)^^(Cgri - 0> ? > 1> / A 1 1 \ 



where Rrad is the radiation amplitude, 

kr sin (5a - Ss)' 

Srad is the radiation phase, 



(A.12) 



Srad = argiF{kr)) + 5s + Sa-^, (A.13) 

and H{x) is the step function, i.e., H = 1 when x >0, and H = otherwise. 

We show that these results are consistent with Eqs. (3.8) and (3.14). Indeed, computing the 
singular contribution from the integral (A.5) in the region ^ <C — 1, we find by similar technique 
that: 

£ W^/'' = (''^'^^^"^' - nK)e-'^') . (A.14) 

Matching this boundary condition with Eq. (2.18) in the region ^ ^ —1, we find: 

\F{kr)\ = —R{5), aig{F{kr)) =Tr-d. (A.15) 

TT 

Let us specify the inhomogeneous solution dU{^;cES,S)/dc at the minimal tail phase S = Ss + 7r/2, 
then the homogeneous eigenfunction ipiiO is excluded from (2.16) and (A.5), since 71 = 0, see 
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below Eq. (2.19). Setting the value 6 = 5s + Tr/2 in Eq. (A.15), we finally find from Eqs. (A.12) 
and (A. 13) that 

Rrad = 2 ■ \^ 7^ = 2R{6a), ^rad = ^a, 

sm{da - 6s) 
where we have used the relation (2.19). 



Appendix B 



In this appendix, we use a local momentum conservation law to derive the dynamical equation for 
velocity c. The local momentum conservation law has the form, 



[1 21 


+ 




.2\ 


t 



1 



-u^ + uuxxxx - UxUxxx + 7:Uxx + o"o'" + "lU u^x + -jasu 



0. (B.l) 



Integrating Eq. (B.l) over |^| ~ 0(1) and substituting the perturbation expansion (3.3) up to order 
O(e^), we derive the following equation. 



UESUid^ ) ^ + ^1 



= 0. 

(B.2) 

The stationary solution uis{^) is given by Eqs. (3.12) and (3.15), and the oscillatory tail in front of 
the embedded soliton is given by Eqs. (3.8) and (3.14). When those formulas are utilized, we get 



/oo 
UEsuid^ = E, 
-00 



(B.3) 



and 



?^+oo 



2kX2kf - l)R\5a) = -r. 



(B.4) 

where E and T are defined in Eqs. (3.18) and (3.19). When relations (B.3) and (B.4) arc substituted 
into Eq. (B.2) and eci(r) is replaced by c(t)— ces (3.4), the dynamical equation (3.20) is reproduced. 
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Figure 1: Tail amplitude r{c,5a) of symmetric nonlocal waves as a function of velocity c. The 
nonlocal waves satisfy Eq. (2.6) with system parameters (4.1). The radiation phase S = Sa is the 
phase of the anti-symmetric linear mode V'2 given by Eq. (2.15): Sa ~ 0.5737. 
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Figure 2: Evolution of the embedded soliton under the momentum-enhancing perturbation (4.4) 
with h = 1.05: (a) numerical solution profile at t = 20; (b) velocity c of the perturbed embedded- 
soliton versus time t, the dashed line is the embedded-soliton velocity ces = 2.64; (c) amplitude of 
the perturbed embedded-soliton versus time t, the dashed line is the embedded-soliton amplitude 
0.9; (d) radiation tail amplitude versus time t. The solid lines in (b,c,d) are theoretical predictions, 
and circles are numerical values. 
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Figure 3: Evolution of the embedded soliton under the momentum-reducing perturbation (4.4) 
with h = 0.95: (a) numerical solution profile at t = 50; (b) velocity c of the perturbed embedded- 
soliton versus time t; (c) amplitude of the perturbed embedded-soliton versus time t; (d) radiation 
tail amplitude versus time t. The solid lines in (b,c,d) are theoretical predictions, and circles are 
numerical values. 
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